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We present a theoretical framework which describes multiply-charged atomic 
ions, their stability within super-intense laser fields, also lay corrections to 
the systems due to relativistic effects. Dimensional scaling calculations 
with relativistic corrections for systems: H, H~, H 2 ~, He, He~, He 2- , He 3- 
within super-intense laser fields were completed. Also completed were three- 
dimensional self consistent field calculations to verify the dimensionally scaled 
quantities. With the aforementioned methods the system's ability to stably 
bind 'additional' electrons through the development of multiple isolated re- 
gions of high potential energy leading to nodes of high electron density is 
shown. These nodes are spaced far enough from each other to minimized the 
electronic repulsion of the electrons, while still providing adequate enough 
attraction so as to bind the excess elections into orbitals. We have found that 
even with relativistic considerations these species are stably bound within 
the field. It was also found that performing the dimensional scaling calcu- 
lations for systems within the confines of laser fields to be a much simpler 
and more cost-effective method than the supporting D=3 SCF method. The 
dimensional scaling method is general and can be extended to include rel- 
ativistic corrections to describe the stability of simple molecular systems in 
super- intense laser fields. 
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I. INTRODUCTION 



The generation of stable, multiply-charged atomic ions via exposure to super- 
intense laser fields is a topic which challenges preconceived notions for ionic atoms 
and is, therefore, of fundamental importance in atomic and molecular physics^. 
Over the past decades, advancements in spectroscopic methods have yielded verifi- 
cation of mono-charged Calcium and Strontium atomic anions^lEl and various gas- 
phase poly-charged molecular ions ^ 7 ^ . However without the large charge volume 
which is provided by the heavy atoms -above- or small molecules it is unlikely that 
species would be able to bind more than one excess electron; this can be noted by 
the relative stability of O -2 in the liquid-phase, yet it's instability within the gas- 
phase^. Theoretical works have developed an absolute upper-limit to the number of 
electrons which may be bound to a atomic center^ N c > 2Z, with N c being the 
number of electrons and Z being the Coulomb charge of the nucleus. Within the 
context of Lieb's frameworks, Hydrogen would therefore be disallowed any excess 
electrons beyond that which yields the Hydride state, thusly H 2 ~ is unstable ^ 1 *! 
Supporting theoretical works have come -some including implementation of 

finite- sized scaling^- and have conclusively determined at gas-phase, dianionic atoms 
are unstable. 

It has been shown that stable, multiply-charged atomic ions may be developed 
within extremely strong laser fields on the order of 10 16 W/cm 2 and above^El. 
Within the field, the electron density - still being bound to the nucleus - has been 
found to be nodal in nature as the Coulomb potential splinters under the influence 
of the field into distinct, localized regions whose positions are governed by the field 
parameters of the laser. This phenomenon is most easily - and best- discussed within 
the context of the Krammers Henneberger (KH) reference frame, electron centric 
frame, where the electron is treated as the stationary body and the nucleus traverses 
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the path of the applied field; in this context the local nodes of electron density are 
located at the turning points on the path of the nucleus. These are the location 
at which the angular velocity of the nucleus decreases and thus spends more time 
in a local area - thus generating a greater pull in that area. Within these nodal 
regions, the bound electrons maintain a great enough distance from one and another 
to minimize their Coulomb repulsion while also giving each a center with which to 
bind. In this field, the electrons -which intuition tells us would be completely ionized- 
are capable of stably binding into multiply-charged atomic ions. The field strength 
allows one to manipulate the location and pull of the nodal centers , thus generating 
a method of control over the potential and therefore establishing the ability to push 
the electrons into and past their most stable state by means of manipulating laser 
parameters, frequency and intensity. 

The contained theoretical works are concerned with High-Frequency Floqeut The- 
ory (HFFT) which allows for a time-independent treatment of the coupling of the 
static Coulomb potential with a time-varying electromagnetic field. This is possi- 
ble by exploiting highly oscillatory fields in which the electrons would be prohibited 
from coupling with the periodic nature of the field due to extremely short periods 
(large frequencies) of oscillation, thus the system's electrons would feel a period av- 
erage of the applied potential coupled with the static Coulomb potential, again this 
is best discussed within the KH frame. This time-average allows generation of the 
aforementioned nodal structure and therefore permits the stability of the subsequent 
states and allows the system to forgo autoionization. The above discussed method- 
ology was introduced to atomic systems by Pont et aP^, van Duijin et al 1 ^ 7 , and 
was used again by Wei et apMSDIHH ^ describe non-relativistic, multiply-charged 
atomic ions. Herein we shall propose a framework utilizing HFFT as a backbone for 
applying relativistic corrections to atomic ions in a time-independent manner. 
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II. RELATIVISTIC CORRECTIONS 



A. Non-Relativistic Methodology 

Consideration within the non-relativisitic cases lies no longer with both the mass 
and the magnetic coupling, but with the time dependent electric field coupling with 
the system's Coulomb potential; this work was performed by Wei et apMEHHH 
and produced stably bound multiply-charged ions for small atomic centers utilizing 
the field parameter (a ), discussed later, and finding detachment energies on orders 
of O.leV to 1.0 eV. The enclosed works, here, expound upon this by adding the 
necessary relativistic corrections to the previous framework. A free electron within 
an oscillating electric field shall undergo oscillatory motions which are governed by a 
coupling to the field; the electron is said to be 'quivering' with a motion defined by 
a trajectory, c?o(£), and a quiver amplitude, «o- A bound electron within the same 
situation shall feel a new potential which is a stacking of the applied field and the 
Coulomb potential of the central charge; the total potential for the system is said to 
be a Coulomb potential dressed by the laser, denoted as a dressed potential, Vdres- 

Under the auspices of HFFT, introduced above and here^, by applying a highly 
oscillatory laser field with an extremely short period the electrons will lack the ability 
to oscillate synchronously with the the applied field. In this manner, the potential 
felt by the electrons is a period average of the oscillatory field, this new potential is 
a dressed potential under the HFFT approximation, V^ s . In all cases addressed 
within this paper the laser-coordinates are: laser fired in y-direction, electronic com- 
ponent linearly polarized in the z-direction, and the magnetic component in the 
x-direction 

The situation of the dressed potential, V£;f s FT , is a time independent problem as 
the field has been period averaged, due to this the full Hamiltonian can be treated 
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within the Time Independent Schrodinger Equation (TISE): 

M = = ^V 2 ^ + VjZTVi. (1) 

Accurate solutions to the equation are difficult for systems of more than one electron 
due to many-body interactions, but approximate solutions can be obtained in a self- 
consistent method (SCF) through Hartree-Fock (RHF/UHF), Density Functional 
Theory (DFT), or post-Hartree-Fock Methods. 

B. The Relativistic Mass Gauge 

Mass -as a fundamental- is conceptualized in two different manners within physics, 
these being the rest mass and the relativistic mass. Rest mass, or invariant mass, 
is for a specific body a constant, measurable quantity denoted m . In opposition to 
the rest mass is the variant quantity: relativistic mass, m r . The relativistic mass 
depends upon the velocity of the observer. The variant nature of m r is a correction to 
the rest mass which accounts for a non-zero kinetic energy for the measured system. 
This means that the relativistic mass increases in magnitude as the velocity of the 
system increases, and shall reach infinite mass as the system reaches the speed of 
light. 

This portrayal of the mass shall be implemented within the Time Independent 
Schrodinger Equation (TISE) for the enclosed work, we shall now need to express 
alterations to the rest mass in terms of the system's laser parameters. 

We are now called to introduce the concept of ponderomotive energy, U p ; this 
being the cycle average kinetic energy of a quivering electron, i.e. electron undergo- 
ing oscillatory motion due to an external field and also qualifying under the dipole 
approximation. This is quantity discussed in context of such systems by Joachain, 
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Dorr and Klystra^ 1 . 



4m P u; 2 



In the above, S Q is the peak strength of the electric field, u is the angular frequency 
of the applied field and both e and m e retain their conventional meanings: of mag- 
nitude of electron charge and electron mass, respectively. This quantity aides in the 
evolution of the rest mass to the relativistic mass as: 

m r = pressed = + 2 _ P N§ /g^ 

m e cr 

= m e (l + 2g)i (3b) 

As can be seen, the quantity q begins to shift the mass and becomes the dominant 
factor within the expression as it approaches unity^. The form of Eq. (3b) was 



found by Brown and Kibble^ and later verified by Eberly and Sleeper^ via the 
Hamilton- Jacobi equation. 

Below we shall discuss the above mass concept as it applies to the TISE for one 
electron (which can then be generalized to N electrons) , these lines shall be discussed 
stepwise. 

Eq.{\) V 2 ^ + V% (4) 

2m e (l + 2q)i 

-> ; ~ k \ 2 2n1 V 2 ^ + VVt (5) 

2m e (l + 2^J#), 

-h 2 

-»• rV 2 *^, (6) 

2(l+2.66xl0- 5 a 2 u; 2 )2 

The first line, Eq. Q, shows the form of the TISE as it appears accounting for the 
mass gauge, which is tuned by the quantity q. Secondly, we have introduced and 
employed the field coefficient, a = ^, as a means of defining q in terms of known 
laser parameters. Lastly, we express all quantities in h=m. e =l units (atomic units), 
this allows us to maintain the relativistic alterations as a unitless multiplier, Eq. (|6]). 



As the multiplier which transforms invariant mass to relativistic mass is a unitless 
quantity, the resultant energies from the final line, Eq. (J6|, shall be in Hartree E#, as 
they would if one ignored the mass gauge entirely. In all cases considered within this 
paper the potential function, V, shall be dressed under HFFT , making y=Vj^f s FT 



C. Trajectory Corrections 



High-Frequency Floquet Theory (HFFT) was first introduced to similar systems 
by Pont et aP^, relies on the frequency of the externally applied electromagnetic field 
to be so quickly oscillating so as even the electrons are incapable of coupling their 
motions to the field. In this manner the D-dimensional dressed Coulomb potential 
-which is in essence a time-dependent problem is simplified to a time-independent 
problem: 

/ \ 



V. 



HFFT_ 

dres 




(7) 



To apply the above period average to a system one must develop an interest in the 
trajectory, a(t), of the laser's path as it's components are required in the above aver- 
aging as ai along with the \ th Cartesian component, X{. Earlier workg^ 15 ^ 18 ^ 19 ! have 
concentrated upon non-relativistic systems, and thus the laser trajectory is equiva- 
lent to the path taken by a free electron undergoing influence by a time-dependent 
external electric field (or laser field where no magneto-coupling is considered); the 
up and down oscillatory motion of the electric field governs the trajectory of the 
electron, seen in Eq.((8j), where the polarization is as discussed in §A. 



a(t) = (a xj a yj a z ) 

= (0, 0, a cos(cut)) . 



(8a) 
(8b) 



Now concerning ourselves with the electronic-magnetic coupling within relativistic 
regimes; this shall be described analogously to the non-relativistic case above by the 
path taken by a free electron in an electromagnetic field now with the electronic- 
magneto coupling accounted. Within a plane-wave laser field, the electron classical 
trajectory can be obtained analytically. For a linearly polarized laser field, the non- 
zero field components can be written as: 

E z = E cosr), (9) 
cB x = E cosr], (10) 

where r\ = out — ky is the phase of the field. By using Newton-Lorentz equation, 

dP 

-^ = -e(E + vxB), (11) 



we have, 



dP _ dP x dP y dP z 

~dt ~ ^~dT , ~dT , ~dT^ (12) 
= {-eE /3 x cosr), 0, -eE (l - f3 z )cosr]} 



where 7 = y 1 + P 2 /(m^c 2 ) is the Lorentz factor. For an electron initially at rest, 



from Eqs.(12 13), we obtain: 

P { Px 1 Py 1 Pz / 



2m e cu 2 u> 
From the above an electron's trajectory can be acquired 



r e*El . eE . ( 14 ) 
{- -stn 7], U, sinrj] 



C 2 . _ x „ C 



a = {-—Q sin(2r}),0, -Q cosr]}. (15) 

oUJ UJ 
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Where Qo = eE /(m e uc) and electron drift motion has been neglected. Within 
dipole approximation, 



which means: 



a = {0, — —Qlsin(2ut) , —Q cos(ut)}, (16) 

oU) 00 



d 2 a , , 

dt (17) 
= {-cw(5o'S^(2wt), 0, — aoQocos(ojt)}. 



Eq. (17) tells us that we can approximately take the above relativistic trajectory as 

an equivalent one for the electron moving in the following effective electric fields, 

-> r um e c ^ 9 , . um e c „ , „ . , 

E={—^Q 2 sin(2ujt), 0, -Q cos(ut)}, (18) 

which will be used in HK theory. In atomic units, the trajectory can be written as, 

a = {—alafsin(2ojt),0,a cos(ujt)} (19) 

By comparing with the non-relativistic trajectory, we have an extra oscillation motion 
along the laser propagation direction, which comes from the magnetic coupling, also 
intorduced is the fine structure constant which mediates the magnetic-electronic 
coupling term, otf. This results in the famous figure-8 motion. To illustrate the effects 
of this correction in HK effective potential, we will first take a 1-D box potential as 
an example in the following. 



D. 1-D Particle in a Box 



For simplicity, we take the 1-D box potential as follows, 



, 7T, Izl < 1, 

V{z) ={ ~ (20) 

0, elsewhere. 
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Then the HK effective potential can be acquired analytically by the following inte- 
gration. 



Veff(z) 



1 

2^ 



2tt 



V(z + a cos(tt))dtt. 



(21) 



When < a < \, 



arcos 
arcos 



eff\ 



z+1 
ao 

z-1 

OLQ 



, -a < z < -1 + a , 
1— a <z<l + ao, 



it, — 1 — a < z < 1 — a , 
0, elsewhere. 



(22) 



When a > \ 



2 ■ 



V, 



eff\ 



arcos 
arcos 
arcos 



z+l 
ao 

z-1 
a 

z-1 



0, elsewhere 



, -«o < z < 1 - a , 
, 1 - a < z < -1 + a , 
, -1 + ct < z < 1 + a , 



(23) 



FigjT] shows how the box potential is modified by the external laser field. It is 
clear to see that as the laser intensity increases above certain value, the original 
potential will evolve into a double well. Moreover, the two wells will become more 
separate and more shallow if ao is further increased. This indicates two important 
features for the ground state in this potential. Firstly, the electrons will become less 
bound or the potential might have higher ground state energy. Secondly, if we have 
two electrons in this effective potential, they will have more space in which to avoid 
each other, which means the electron repulsion energy will tend to be smaller for 
un-paired electrons. Hence the final ground state energy for multi-electron systems 
should depend upon the competition of these two factors. 

Once we introduce the relativistic corrections to the electron trajectory, the effec- 
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tive potential along z-axis will becomes: 

V efJ (z)=^ Jd{l(^{z + a costt) 2 +a 2 a 2 sin 2 2tt'j , (24) 

in which the integral will be calculated numerically. From Fig. ([2]), it is interesting 
to note that, when a is large enough, for example a = 10 as in Fig|2^b), the 
effective potential will have three local minima. Qualitatively, this should result 
from the relativistic figure-8 motions, which helps the electron maintain a position 
nearer the orbit center for a greater period of time. Another characteristic is that 
the three separate minima become much shallower for higher laser intensity. These 
drastic changes over the effective potential will be expected to have influences over 
the many-electron states bound by the potential. It seems that the relativistic effect 
provides us another way to engineer the potential. Based upon the observations 



over Eq.(24), we can even think about using lasers with different colors superposed 
together to have more freedoms modifying the effective potentials. Fig. [3] shows us 
one example with two alternating electric fields as follows, 

a{t) = a cos(ut)e x + aisin(4:Ut)e z . (25) 



E. Potential under Consideration 

Beginning with the spherically symmetric Coulomb potential, the applied external 
fields shall morphologically alter this potential to conform to the trajectory discussed 
in § B. In doing this there are three main regimes in which the potential may exist: 
firstly, the spherical; secondly, the pseudo-linear; and finally, the parametric. The 
potential only exists within the spherical regime when there is no applied external 
field, as the field begins to evolve the coupling of the external field and the Coulomb 
potential becomes apparent with the electric field component dominating; this creates 
a regime where the system maintains an almost linear behavior as if there were no 
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relativistic correction to the trajectory. Distortion of the pure linear nature exists 
but is a small effect compared to the primary electric effect. 

As the intensity of applied field increases the magnitude of the magnetic field 
begins to compensate for the dampening effect of the fine structure constant, a/; as 
this takes place, the magnetic contributions to the field coupling begin to dominate 
the system, whose character now exhibits the hourglass figure of the parametric 
regime. Figure [4] shows a series of contour plots of the potential energy plotted 
in the x-z directions for a series of field intensities, ao; the behavior of the system 
and it's development through the previously discussed regimes should be obvious. A 
three-dimensional plot of the potential energy surface for o.q = 100 is also enclosed as 
Fig. |5j this plot only considers displacements in the spatial x-z directions for reasons 
introduced in § A, where the vertical axis describes the magnitude of the potential 
energy at this x-z location. 

This potential is ideal for attempting to develop mutliply-charged ionic system 
from small nuclei, as it maintains a potential well at the center of the systems along 
with four other locations along the parametric trajectory; this allows bound electrons 
to attach to the individual potential wells while maintaining a large enough distance 
from each other to minimize electron-electron repulsions. 

III. METHODOLOGY FOR D=3 CALCULATIONS 

Three dimensional calculations which describe our systems were executed as a 
mean of verifying the simpler Dimensional Scaling approach discussed later. The 
methodology consisted of unrestricted Hartree-Fock (UHF) utilizing the Pople- 
Nesbot equations -which allow for the accommodation of basis sets- to complete 
calculations for a series of total electron counts, N, per single value of the field inten- 
sity parameter, a . The calculations were dependent both upon appropriate choice 
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of basis set and upon the locations selected as the centers for these basis sets; for 
the determinations of the locations of probable electron density we deferred both to 
'exact' locations of the electrons from the limit D— > oo calculations and to contour 
plots of the potentials for a given field intensity value, see Figure [4] for an example. 
Upon discerning from the above information the locations of the psuedo-centers 
within the system-space, a basis set was selected which could describe these nodes. 
There exist, at minimum, 5 distinct nodes of electron density within the system, 
these being at the center coexisting with the origin of our coordinates systems- and at 
four psuedo-centers residing upon the parametric curve described by the relativistic 
trajectory used within the HFFT Potential. 
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X 




< y 


> = < 


z 









-(*)«n(2t) 


aocos(t) 



(26) 



The orbital centers were selected to satisfy the above curve and to coincide with 
the "hottest" locations displayed within the contour plots of potential energy. A 
basis set was selected which allowed for significant description of both polarized and 
diffuse phenomenon residing on small centers. The centers of potential electron den- 
sity which do not reside at origin were described by basis sets optimized to describe 
the atom of the nucleus at the system's origin, i.e. all five centers on hydrogen use 
a hydrogen-fit basis set. Overlap and kinetic integrals were performed with varia- 
tions on the integrals described by McMurchie and Davidson in their seminal paper 
here®'. Numerical integration methods were appropriate for the more challenging, 
non-analytic potential integrals. A global adaptive method was applied to the sys- 
tem which performs admirably with Gaussian functions placed at the coordinate 
origin; but as the method relies on sampling points within the equation-space to find 
non-zero areas of the function, Gaussians placed distances away from the origin were 
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sometimes so small relative to the distance from origin as to be invisible. For this 
reason two-centered intregral includes the cost of shifting the coordinate-origin to 
the center of the product Gaussian as defined by^: 

4>ff(a A , r-R A )x(pg F (a B , r-R B ) 

= exp{-(a A (r-R A ) + a B {r-R B ))} 

(27) 

= exp{-a a a B (a A +a B ) 1 \R A -R B \ 2 }x 
<t>? s F {at A +a B ,r - R P ). 

In this way the chosen method of numerical integration was capable of adequately de- 
scribing the three-dimensional potential energy integrals. Single electron cases were 
verified prior to enacting the self-consistentt field calculations, as the energies of the 
single electron system may be revealed as the eigenvalues of the canonically orthog- 
onalized H core matrix alone, H^ e =T MjJ ,+V Ati j,. The four-centered integrals needed 
for the self-consistency were generated by exploiting the axillary function defined 
by Boyes 26 . The calculations for the mutli-electron energies were performed self- 
consistently with a convergence set to six decimals of accuracy, as chemical accuracy 
is define as 1.6mE#, this set limit should suffice. 

A plot of the square of the linear combination of atomic orbitals which comprised 
the set describing the appropriate eigenvalue yields semblance to the wave function of 
the system, whose probability density(|\l/| 2 ) is shown in Figure [7] for H~ and for He - , 
a two and a three electron case. It should be noted, as the D=3, UHF calculations 
were employed to verify the results of the Dimensional Scaling approach, that the 
locations of electron density shown in the aforementioned probability density plots 
speak to the validity of the Dimensional Scaling approach as the locations of the delta 
functions, see discussion provided in Results Section, which describe the seemingly 
stationary locations of the electrons in the limit of D — > oo calculations are very 
near those locations of highest probability density for finding elections given by the 
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D=3 calculations. Three-dimensional verification was performed for both the H 
and He~ species, as the regions of high potential energy become very delocalized 
for the remaining high-field species. Aside from the high level of derealization 
prohibiting the species from being described sufficiently with an appropriate number 
of Gaussians, the potential also spans a region of space on the order of 100 Bohr radii, 
yet optimized Gaussian basis sets for atomic centers span an order of 1-10 Bohr radii; 
this prohibits overlap for these species forcing the use of more and more Gaussians. 
This acts prohibitively as the matricies required for UHF calculations are nxn, where 
n scales as the number of centers by the number of basis contractions; this obviously 
limits the achievable intensities applied to systems which can be calculated in this 
manner, especially on stand alone machines. 

IV. DIMENSIONAL SCALING: CALCULATIONS AND 
CONSIDERATIONS 

A. Methodology 

Many Body interactions are something which has troubled computational method- 
ologies within quantum mechanics since inception; throughout the years the physical 
and chemical communities have made great advances in the field of electronic struc- 
ture theory to help account for these electron-electron interaction through variational 
practices such as the Hartree Fock Method or Density Functional Theory. The alter- 
native method to the aforementioned is a Dimensional Scaling treatment pioneered 
by Herschbachp3, discussed irP&E2 anc [ j s briefly introduced here for the central force 
problempS. Given the TISE for the simple central force problem in D-dimensions: 

[- l -V 2 D + V{r)]^ D = e D ^ D . (28) 
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If we were to employ polar coordinates to the above we would require: 

D 



k=l 



1 d 



D 



D 



dr 



dr' 



T 2 



(29a) 
(29b) 



Where Eq.(29a) gives the definition of the radial coordinate in a generic D-scaled 
space, and Eq.(29b) is the polar Laplacian in this D-scaled space, h 2 D is the term 
which retains all angular dependencies. These angular and radial terms shall be dealt 
with in an divide and conquer treatment reminiscent to the radial and angular terms 
of the Rigid Rotor /Harmonic Oscillator approximations for the simple diatomic. We 
first write the wave function in D-dimensions to be the product: ^ d = rC y{^D-i), 
where all radial dependencies are in the r' term and the D-l remaining angular 
dimensions are described through y(Qo-i)- Now solving the angular terms for the 



form Eq.(30), and the recognizing that the V(r) term in Eq.(28) can be set to equal 



shown in Eq.(31 ). 



magnitude as the ejj term, thus making Eq.(|28j) reduce to the Laplace equation 

L 2 D _ x y(si D - x ) = cy{n D ^) (30) 



V>'^D-l) = 

W + D-2)-C}r t - 2 y(Q D _ 1 ) 



(31a) 
(31b) 



This means: C=C(( + D — 2); and the Hamiltonian Operator in Eq.(28) is now of 
the form: 

C{£ + d-2) 



H 



D 



-K D -i(r) + 



+ V(r). 



(32) 



2 "" *v > ■ 2r 2 

In the above, Ko_i(r) is the single non-angular term from the polar Laplacian in 



Eq.(29b). We may now pass the system through a unit Jacobian, making: JdI^dI 
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$£0 Jd is the radial part of the unit Jacobian and is of the form: r^ D x ). This would 
mean, = r~2^ D ~ 1 )^ D , All leading to the form of Kp_i seen here: 

*WO=r * ( H-^— Y~Y-7T} ( 33 ) 



Reassembling all the above, and placing them appropriately back into Eq.(28), 
one shall -after menial simplification- get: 

{-^ + ^ + ^>>*-^ (34, 



Eq.(34) is the radial, D-scaled form of Eq. (28), where the only dimensional de- 
pendencies lay within the A terms as: A = t + \{D — 3). The above leads to the 
minimization problem defined by the Hamiltonian discussed in §B. 

B. Planar Infinite-D Hamiltonian 

Prior works publishe d 19 1 114 1 have also described systems both by infinite dimen- 
sional limit and then verified with three dimensional self consistent methods. The 
dimensional scaled Hamiltonian presented in previous works was diatomic in nature 
and of the formpSl 

N , N 



N N-l t (35) 

+ E E / 

<=1 j=i+iy i z i- z j) 2 +Pi+P 2 j 

and relied on previous works in which D-scaled Hamiltonians for diatomic systems 
were constucted, these diatomic Hamiltonians are also of the form above, and denoted 
Hda- 

Hamiltonians of this form are applicable to the previous works as those non- 
relativistic systems, this is due to the consideration in absence of the second degree 
of symmetry breaking in the linear potential systems. 



The diatomic-based Hamiltonian performed well on linear systems, but when at- 
tempting to use the above described Hamiltonian on a relativistically corrected tra- 
jectory, it was found that the equation behaved erratically with a smoothly evolving 
intensity as seen in Figure [6j and was in disagreement with the early 3-D calcula- 
tions. Although this Hamiltonian does not work well, overall it was found to be in 
good agreement while the potential was in the spherical and pseudo-linear regimes, 
failing only as the system approached and entered the parametric regime. 

As dimensionally scaled Hamiltonians are not unique in nature, they are not 
singular in form. We relied on arguments based on the breaking of spherical and then 
cylindrical symmetries to generate the followed assumed Hamiltonian for systems 
with broken radial and cylindrical symmetries, yet maintaining three orthogonal 
planes of symmetry. 

As the dimensional-scaled, single electron central force problem yields as it's 
Hamiltonian Wcf, 

KcF = ^r 2 +V(r). (36) 

One can see from above the spherical nature of all terms within the Hamiltonian, 
as the potential is radial. This equation predicts the ground state Hydrogen energy 
to be at —^Eh, exactly where it should be, and predicts the inter-atomic distance 
between the electron and the proton to be 1 in units of Ksohr- If one were to alter 
the above potential to conform with either the relativistic or non-relativistic cases 
discussed in this paper, the energies obtained would possess no physical significance 
and would overall behave similarly to the diatomic equation with the relativistic 
trajectory. The addition of multiple electrons to the spherical symmetric problem 
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yields : 



N N 

n Z\ 1 (37) 

+ E E -r=^- 



Now we examine the diatomic case shown in Eq. (35). This equation has been 
designed to allow for potential which are of a cylindrical nature, as it must for a 
diatomic system. Reevaluating the equation in Cartesian coordinates yields: 

1 N 1 N 

Z i=l X i + % i=l 

N 2V-1 j (38) 



Reevaluation of the spherically symmetrical case in Cartesian coordinates gives: 

1 N 1 * 

^ = 2 E x 2 + 2 + z ? + E ^ y*> z ^ 

N N-l (39) 

+ E E / • 

t=i 3=i+i \JA + x 2 j + yl + y] + ^ + z] 

As can been seen from the above, the orthogonal coordinate which is being cleaved 
from spherical symmetry has been removed from the kinetic evaluation and have been 
treated as a difference in the electron-electron term. Continuing to use Cartesian 
coordinates, as it is this coordinate system which makes the relations apparent, we 
can move to an equation where the x-coordinate is now allowed to deviate from 
radial symmetry. This again would remove the symmetry breaking coordinate from 
the kinetic term and utilize it as a difference in the electron-electron term. This 



Hamiltonian is shown in Eq. (40). It's energies completely agree with those of 
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the the radial and cylindrical cases, and by allowing this symmetry breaking in the 
x-coordinate can be use for the potential discussed in Section II §E. 
This Planar Hamiltonian, Hp: 

N „ N 



Up =\ ~2 + v ^ Vi > z ^ 



2 ^ v 1 

»=i y * i=i 

N N-X , (40) 



i=l j=i+l J (z-Zj^+ix-Xj) 2 -^^] 

was found to allow -but not require- the breaking of symmetry, as it was applied to 
the previously discussed linear systems and had extremely good agreement, shown 
in Figure [6] The relation between the Cartesian coordinates used in the above 
and the geometry of the system is shown in Figure [8] When applied to our system, 
containing the corrected trajectory, this new Hamiltonian performs both smoothing 



with smoothly evolving trajectories- as seen in Figure [10] This planar Dimensional 
Scaled Hamiltonian was capable of adequate description of the system in all three of 
the regimes discussed earlier. 

V. RESULTS AND DISCUSSION 

We shall concern ourselves with a discussion of the binding energies (B.E.) for the 
following species: H~, H~ 2 , He~, He -2 , He -3 ; where the B.E. is the negative of the 



detachment energy for a single 'excess' electron shown in Eq. (41), where N signifies 
the number of electrons for a species. 

B.E. = Eh(N) — E H (N — 1) (41) 

Figure [9] Displays the binding energies for the two species, Z=l and Z=2. From this 
figure we can see a clear maximum binding energy for the H~ species (top left) at, 
roughly, a = 10; this energy shows a stability of the second electron of 0.047 Hartree 
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(1.28 eV). Also shown in Figure[9]the B.E. curve for H 2 ~ (top right), showing a stable 
binding of 0.00012 Hartree (0.0033eV). This, the B.E. of H 2 ~, reaches it's minimal 
value asymptotically with increasing a implying the addition of any further electrons 
will not be allowed; this fact was verified by performing the requisite minimization, 
and if the mass gauge was not applied to the system the number of allowed additional 
electrons would increase unrealistically and seemingly without bound as the laser 
intensity is increased. 

The middle, left plot in Figure [9] shows the binding of a third electron to Helium 
at a equals, roughly, 10 again; the binding energy for this species at it's greatest 
magnitude is extrapolated to be 0.057 Hartree (1.55 eV). The second 'additional' 
electron to Helium (middle right) is most stable at «o=90, with a B.E. of 0.007 
Hartree (0.19 eV). The fifth and final electron which can stable bind to a Z=2 center 
(bottom center) is, similarly to the H 2_ species, a terminal binding who reaches an 
asymptotic stability with increasing alpha, the B.E. is 0.0004 Hartree (0.011 eV). 

Of the two relativistic corrections accounted for in the above framework, the 



trajectory is the paramount addition. An examination of Figures [10] and [6] shall 
be required for the subsequent discussions. The introduction of the planar Large-D 
Hamiltonian, Hp, for the systems was a boom which aligns itself in agreement with 



the previous works -see Figure 10 a and compare plots a and c in Figure [6} and yet 



out preforms the previous equation when this relativistic trajectory is introduced, to 
see this compared plots b and c in Figure [6] 

By comparing plots a and d of Figure [6] it can be seen that the raw energies and the 
binding energies between the non-relativistic trajectory (plotted with Tin a) and the 
relativistic trajectory are not extreme. This can be verified and more clearly seen by 



examination of Figure 10 c, wherein both the non-relativistic and relativistic binding 
energies are shown; although the energies become quite different with increasing qjq, 
the field intensity which yields the most stable binding energy is same and the most 
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stable binding energy deviates only in the thousandth of a Hartree. Consideration of 
the mass gauge for this system provides a very slight correction within the values of «o 
examined here; where although these values of «o indicate laser field strengths on the 
order of atomic units and greater (in competition with the Coulomb potential of the 
center) they are in no way strong enough to generate quivering electron masses, m r , 
which deviate significantly from the invariant mass, mo. Due to this, the deviation 
of the binding energies due to consideration of m r over m is also slightly less than 
breath-taking, this can be seen in Figure [10b and again impressed by examination 



of the deviation of the raw energies of Hydrogen in Figure 10 i. To be gleaned from 
this is that within the examined field strengths, m r deviates very little from mo, but 
more significant is the fact that the overall consideration of relativistic effects does 
destabilize the system, but not by an appreciable amount. 

Figure [12] helps emphasis the differences between the corrected and non-corrected 
systems. The top level of the figure is the probability density, \ip\ 2 , of corrected (left) 
and non-corrected (right), the density in the outer orbital centers is lower due to a 
more diffuse spread of probability density at these locations compared to the non- 
corrected. The mid-level of the plot is the negative probability density, superimposed 
on the probability density function is the trajectory of the system which it describes. 
Below all else is a contour plot of each system, again emphasis is merited on the 
more diffuse spread in the corrected system due to the evolution of the parametric 
trajectory as opposed to the linear oscillating trajectory of the non-corrected system. 



Figure 11 shows a plot of the B.E. For the H (left) and He (right) systems 



from both the SCF (blue) and D-Scaled (red) methods. The lower B.E. plots in 



Fig. [TT] have been normalized to the minimum value to show the tight agree- 
ment between the qualitative assessment of the scaling procedure and the SCF 
method. The Dimensionally Scaled minimization problem bore 'exact' position of 
the electrons as (x u y x , z X ] a )=(4.1660 x 1(T 10 , 5.4461, -12.2387; 20) for the Z=l, 
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N=l system and (x u yi, z x \ x 2 , y 2 , z 2 ] a ) = (6.4472 x 10" 8 , 5.0386, 16.5852; -6.4472 x 
10~ 8 , 5.0386, —16.5852; 20) for the Z=l, N=2 system; these localized electron posi- 
tions are similarly predictable as they attempt to bind to regoins where the angular 
velocity of the nucleus is lowest in the KH Frame; these locations would most notable 
be the |'s the trajectory. With a single electron the central point of the parametric 
curve, set a origin in our calculations, binds the electron strongly; as more electrons 
are introduced they are situated at locations which minimize the electron-electron 
repulsion of the system. A Mulliken population analysis of the systems shows that 
the orbitals about the central charge typically possess a smaller number of electrons, 
except in the N=3 (and assumably the N=5 case); the results of such population 
assessments can be seen in Table [TJ In this way we are able to verify not only the 
energetic behavior of the D-Scaled Hamiltonian but its treatment of the electrons in 
space. 

VI. CONCLUSION 

It has been shown that even under conditions of relativity multiply-charge atomic 
ions should still be achievable within the confines of an intense laser field. The sta- 
bility of several atomic-anions has been found and discussed, some ions are on the 
order of an entire electron volt more stable than the ionized system. The importance 
of the general dimensional scaling procedure was verified not only via the energetics, 
but with a comparison of the 'exact' locations of the electrons as predicted by Dimen- 
sional Scaling as they compare to the probability densities from the standard SCF 
procedure. These species were found to be stable and should, therefore, be experi- 
mentally realizable. Stability of simple molecular systems in super-intense laser fields 
have been previous discussed here^HSSl. This dimensionally scaled framework with 
relativistic corrections yields itself easily to a description of molecules and molecular 
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ions within the confines of super- intense laser fields, which shall be undertaken next. 
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VIII. APPENDIX 



A. Elliptical and Circular Polarizations 

Within the above concern was only given to the relativistically corrected linear 
polarized light, this is because similar corrections placed on elliptical polarized light 
yield no new -or interesting- phenomenon. Following similar mathematics as to 
achieve an analogous relativisticly considered trajectory would consider the electronic 
coupling as before and consider a minor weighted coupling to the magnetic field, as 
in the linear case before. An elliptically polarized laser fired in the y-direction with 
the electronic major (minor) axis oriented in the x(z)-direction yields the trajectories 

«s,BHi P w(*) = Ui c °s(<P)i 0, e 2 sin{(f))} (42) 
® b, Elliptical) = {-P2Sin(<P), 0, p lC os(4>)}. (43) 



Within Eqs. (42 43) the amplitude in each the major and minor axis is denoted 
by the subscript 1 and 2, respectively, and the e and (3 are the coefficients of the 
electronic and magnetic components. The trajectory generated by both the above 
biases applied to free particle merely generates a new ellipse with a major and minor 
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amplitude mediated between those of the above and tilted by an angle respective to 
the coefficients. 



a = {eicos(t) — /3 2 sm(t), 0, e 2 sm(0) + (3iCos((p)}. (44) 



This can be seen graphically in Fig. 13 which displays the individual electronic and 
magnetic trajectories and then in Fig. [14] which displays the combined trajectory 
which would be followed by a free particle traveling within this electromagnetic field. 

If this scheme is applied to a circularly polarized field the same mathematics will 
appear but the coefficients within the trajectories will be reduced to ei = e 2 = e and 
0i — P2 — P- This will yield the same uninteresting phenomenon, but merely present 
it as the mediate of two circles with no change in the orientation angle as there is no 
unique point of reference on a circle. 
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FIG. 1: Non-relativistic effective potential for a ID particle in a box under different 

laser intensity, measured by a . 
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FIG. 2: Relativistic corrections to the effective potential for different laser fields. 
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FIG. 4: Contour Plots of Dressed Potential for (clockwise) a =0, 25, 100, 50. Note 
both the shift in regime as «o grows and the key below the plots for the 
interpretation of the intensity of the contours. 
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FIG. 5: Three Dimensional Plots of the Potential Energy, V£;f s , as a function of 
x and z coordinates for the case ao=100. Left and Right of above are two different 

angles of same surface. 



Species 


Inner Oribital 


Outer Orbital 


Total Electron 




Population 


Population (each) 


Count (N) 


H 


0.01645 


0.24576 


1 


H 


0.000524 


0.49987 


2 


He 


0.000432 


0.499892 


2 


He- 


1.00001 


0.499999 


3 



TABLE I: Results of Mulliken Population Analysis, note that there are four outer 
orbitals the table contains one of the four values. 
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FIG. 6: Plots of both moecular energies and binding energy for (clockwise): T-Lda 
with non-relativistic trajectory; Tin A, relativistic trajectory; Hp, relativistic 
trajectory; and Hp with non-relativistic trajectory All may be read as 
Yellow: Binding Energy; Blue:Hydrogen Energy; Purple:H~ Energy 
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FIG. 8: The above displays the relationship between the system's geometry with 



respect to the electrons and the coordinates use in Eq. (40) 
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FIG. 10: Plots of binding energy comparisons for (clockwise): comparison plot of 
H.DA versus Hp, non-relativistic trajectory; Hp with relativistic trajectory m and 
m r , differences value of Hydrogen energy between use of m and m r ; and a 
comparison of relativistic trajectory (both mo and m r ) with non-relativistic 

trajectory both using Hp. 
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FIG. 11: Plots of binding energy for the H~ (left) and He~ (right) systems. Top, 

Non-Normalized plots of the calculation data showing agreement between the 
methods. Bottom, -B.E./B.E. maa , to emphasis the qualitative simularity between 
the methods as they share minima for the B.E. curves. 
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FIG. 12: Top Row: Plots of the Probabiltiy Distribution for the corrected (left) 
and non-corrected (right) H~ system, directly below is a superimposition of the 
trajectory upon the the probabilty density function plot to emphasis their relation. 

Bottom Row: Contour Plots of the H~ system, both corrected (left) and 
non-corrected (right), note the different scales on the vertical (z) axis and the more 

diffuse behavior of the corrected system. 
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FIG. 13: The electronic elliptical contribution can be seen as the function whose 
major axis lies in the z-direction (vertical) and the magnetic contribution has an 

orthogonal orientation. 



41 





FIG. 14: Within the above plot, the total trajectory can be seem, the amplitude of 
the major and minor axises are mediated in value between those from the electronic 
and magnetic components and the orientation is set off by an angle whose value 
respects the same coefficients as the relative amplitudes. 
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